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OO . Abstract: We study the profile of a tliin disk configuration as described by an 

axisymmetric ideal magnetohydrodynamics steady equilibrium. We consider the 

["t"] . disk like a differentially rotating system dominated by the Keplerian term, but 

^jH ' allowing for a non-zero radial and vertical matter flux. As a result, the steady 

state allows for the existence of local peaks for the vertical velocity of the plasma 
particles, though it prevents the radial matter accretion rate. This ideal mag- 
netohydrodynamics scheme is therefore unable to solve the angular momentum- 
O ' transport problem, but we suggest that it provides a mechanism for the generation 

.«_> , of matter-jet seeds. 
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I ■ 1. Basic Statements - In this paper, we consider an axisymmetric ideal magnetohydro- 

f~^ I dynamics (MHD) model for the disk structure surrounding a compact astrophysical object. 

■^ ■ Our approach, which differs from the standard astrophysical model for the disk morphology, 

•^^ I is justified in view of the microphysical properties of the plasma having negligible value of the 

t~^ ■ viscosity and the resistivity coefficients. Indeed, the standard model pj [2J accounts for the 

^D I angular momentum transport by postulating a non-zero viscous stress tensor appearing in 

^^ ■ the equilibrium as an effect of the internal turbulence. Despite this model is rather successful 

in explaining the observations, it remains affected by the non trivial shortcoming that the 
microscopic estimations for the viscosity coefficient give much smaller values than those ones 
required to fit data. This is particularly true in accreting objects like X-ray stars, where the 
«2j ' observed rate of accretion is very high. An analogous difficulty concerns the proper balance 

C^ ' of the force acting on electrons, when a significant infalling radial velocity occurs. The solu- 

tion to this puzzle is offered by the presence of an effective resistivity of the plasma, which 
compensate the azimuthal Lorentz force acting on the electrons. Also such values of the 
resistivity coefficient are in contradiction with the microscopical estimations related to the 
astrophysical scales. The turbulent viscosity and resistivity coefficients can be maintained 
both relevant for the equilibrium by postulating a Prandtl number of order unity for the 
accreting plasma. This shortcoming of the standard model is still open to the scientific de- 
bate and leads to speak of an anomalous visco-resistivity of the disk plasma, rising from the 
so-called magnetorotational instability described by E. Velikhov [3,4j- However, the labora- 
tory experience on turbulent plasmas does not confirm the correlation between the internal 
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instabilities and the applicability of a viscoresistive MHD model with Prandtl numbers of 
order unity. 

In our analysis, we address a different scenario as closely related to the transport profile of 
magnetically confined laboratory plasmas [5]. Indeed, in [B], has been shown how the ideal 
MHD profile of a thin disk can have a significant deviation from the Keplerian behavior for a 
strong vertical Lorentz force. In [7l|8] it has been demonstrated that the local configuration of 
a thin disk can be reduced to the profile of a ring sequence, when the electromagnetic backre- 
action is sufficiently intense and the vertical Lorentz force is able to confine the disk (see also 
[3]). Despite a solution of the angular momentum transport {i.e., a non-zero radial accretion 
rate) is not jet available in this ideal scheme, it suggests the idea that (instead of a turbulent 
behavior) the disk plasma could be characterized by the formation of microstructures which 
might drive the typical ballooning-mode instabilities observed in laboratory systems. The 
accretion mechanism would then rely on the porosity of the plasmas near the x-points of the 
magnetic field configuration. 

We here follow the research line traced in [3 [71 |H] , but including vertical and radial matter 
fluxes in the ideal MHD analysis (see also [10l[Tll[T2]). The resulting equilibrium configura- 
tion is characterized by a vanishing net mass accretion rate, as in the previous approaches, 
since contributions due to the radial velocity over and below the equatorial plane cancel 
out. However, our local configuration permits the existence of specific disk layers where the 
vertical velocity can acquire a very large amplitude. This scenario suggests the possibility to 
realize matter jets from the disk equilibrium profile, especially in view of the implementation 
of this same morphology in presence of the expected plasma instabilities and this issue consti- 
tutes the main achievement of our analysis. Indeed, the link between highly magnetized disk 
plasmas and the jet formation (observed in the different electromagnetic bands) constitutes, 
in the standard model, one of the most relevant motivations for studying the angular mo- 
mentum transport across these axisymmetric systems (for a proposal on the jet formation see 
[13]). Our result appears particularly appealing because the viscoresistive model has some 
difficulties to create the conditions for the jet formation, because the viscoresistive MHD 
scenario predicts diffusive magnetic configurations [HI [TS] . 

2. Vertical and Radial Equilibria - To adapt the MHD equations to the axial symmetry, we 
introduce cylindrical coordinates {r, 0, z} and the fundamental variables have no dependence 
on (f). The gravitational potential of the central object writes 



x{r,z)^GMs/ Vr^+z^, (1) 

where Ms is the star mass and G the Newton constant. The continuity equation, in the 
steady state of the disk, reads as the vanishing divergence of the matter flux ev (e being the 
mass density and v the velocity field) and it admits the following solution: 

ev = -er{d,e)/r + e^eruj{r,z^) + e, {drQ)/r , (2) 

where, w(r, z^) is the disk differential angular velocity (in view of the co-rotation theorem 
|16| . it can be written as uj = uj{^p), ip{r,z'^) being the magnetic flux surface). Here, the 
generic function 0(r, z) must be odd in z to deal with a non-zero mass accretion rate of the 
disk (Md): 

/+Ho 
dz evr ^ 47r Q{r, iJ) > , (3) 

-Ho 

where Hq is the half-depth of the disk and Vr < the infalling radial velocity. The magnetic 
field reads as 

B = -er {d,i;)/r + e^ /(^, z)/r + e, {dr^)/r , (4) 

with the function /(V", z) allowing for the existence of poloidal currents in the disk plasma. 



In order to investigate the effects induced on the disk profile by the electromagnetic reaction 
of the plasma, we consider a local model of the equilibrium around a radius value r = tq. 
Thus, we split the mass density e and the thermostatic pressure p into the background (barred 
terms) and perturbation (hatted terms) components, as e = e(ro,z^) + e and, of course, 
p = p(rQ,z'^) + p, respectively. The same way, we express the magnetic surface function in 
the form ip = V'o(''o) + V'i(''Oi i"~i"Oi ^^)i with tpi <C ipo- In particular, the co-rotation theorem 
allows the representation a; = cjk + '^qi/'Ij where Wk is the Keplerian term and ajg = const. 
This form for uj holds locally, as far as (r — ro) remains sufficiently small so that the dominant 
deviation from the Keplerian contribution is due to ipi. 

When we consider the "drift ordering" for the behavior of the gradient amplitude (i.e., 
the first-order gradients of the perturbations are of zeroth-order while the second-order ones 
dominate) we are able to cast the equilibrium configuration by a hierarchical ordering of 
the different contributions. Accordingly, the profile of the disk toroidal currents J^ and the 
azimuthal component of the Lorentz force F^ have the following form: 

J^ ~ -c (a^Vi + 92i^i)/47rro , (5a) 

F^ ~ {dj dri^ ~ dr Id,ij)lA^ro . (5b) 

We now fix the equations governing the vertical and the radial equilibria of the disk by 
separating the fluid components from the electromagnetic backreaction (as discussed in [71 
[8]). Indeed, we assume that the contribution of the toroidal magnetic field does not affect 
these equilibria and therefore the resulting configuration equations overlap that ones derived 
originally by B. Coppi in 2005. The quantity /(^, z) plays a crucial role in the azimuthal 
equation only, allowing for a non-zero radial velocity even in the ideal MHD scenario. 

The splitting of the MHD equations gives, for the vertical force balance, the following 
system of equations 

D = e/eo = exp[-zVi?o] , ^o = AKjiT/{mul), (6a) 

d,p -f ulze - d.iPi {dl4>i -f dl^i) /Anrl = , (6b) 

where eo = eo(''o) = e(''Oi 0) and the temperature T admits the form {K^ being the Boltzmann 
constant) : 

2KBT = m{p/e)^mip + p)/{e + e) = 2KB{T + f) , (7) 

The radial equilibrium splits into the decomposition for w (i.e., oj ~ u)ji + 6u) ~ ujo(Tpo)+'^o4^i) 
and the equation 

2uj^ro{e + e)Cjoi^i + dripi [pl^i -h (9^Vi)/(47rr2) = 

= 9,[p + (9,^i)V(8^r2)] + 9,^1 92^1 /(47rr2) . (8) 

Let us now define the dimensionless functions Y , D and P, in place of i/ii, e and p, i.e., 

Y = ko^i/{dro^o), D = l3i/eQ, P = l3p/po, (9) 

where po = 2KB%eo (To being f{r, 0)) and 

P ^ 8^Po/i?o. = l/(3e') ^ klHl/i . (10) 

Here, we introduced the fundamental wave-number /cq of the radial equilibrium, defined as 
/cq = iuj'^/v\, with v\ = B'^Q/^ATreo), recalling that B^o = (9ro^o)/^o- Thus, we introduce 
the dimensionless radial variable defined as x = kQ{r ~ vq), while the fundamental length 
in the vertical direction is taken as A = ^fe^Ho, leading to define u = zj IS.. This way, the 
relation D = exp[— u^e^] holds. 



These definitions allow us to rewrite the vertical and radial equilibria as follows 

a„2F + e,b - 2 {dlY + e.dlY) d^^Y = , (Ua) 

{D + b/(3)Y + dlY + t^dlY+ 

+ \d^P+{dlY + e,dlY)d^Y^Q. (lib) 

The system above fixes two of the three unknowns y, P and D and has to be completed by 
the azimuthal equilibrium and the electron force balance. We stress that, in setting eas.([TT|). 
we have neglect the inertial contributions coming fi'om the radial and the vertical components 
of the matter velocity, because of their small amplitude. 

3. Additional Configuration Equations - The equilibrium along the toroidal symmetry of 
the disk writes 

evrdr{ujr) + evzdz{u!r) + eujVr — F^f, . (12) 

Using the co-rotation theorem, i.e., lo — Ld{tp), and eqs.(|4]), (j5bl) . the equation above restates 
as 

er^lvrBz — VzBr) + 2eVrUjK/'^o = 

= [dzl {dr„i>0 + dri>l) - drld.^i] / {^^rlCj^) ■ (13) 

In this scheme, we link the radial and vertical velocity fields to the poloidal currents. Indeed, 
we aim to set up a picture in which the radial matter flux is induced by the poloidal currents, 
i.e., establishing a direct relation between Vr and F^. To this end, as well as for the consistence 
of the model, we have to analyze the electron force balance equation: E + (u/c) /\B = 0. The 
balance of the Lorentz force has radial and vertical components which express the electric 
field in the form predicted by the co-rotation theorem [161 IZ] and, since the axial symmetry 
requires E^ = 0, the (/)-component of the this equation retains the form 

VzBr - VrBz = . (14) 

We now observe that, substituting ea.([T^ in ea. (fT5)) . we get the searched relation between 
Vr and i^0, i.e., 

2evruj^ = [dzl (a,„Vo + aW^i) - drIdz^i]/iAnrl) . (15) 

This relation, together with ea. lfTl)) . provides a system for the two unknown functions 9 and 



I, once tpi and e are given by the vertical and the radial equilibria, i. e, eq. (|llap and eq. (|llbp , 
respectively. 

Substituting Vr expressed by ea. ((T5|) into the electron force balance (ITl)) . we get the fol- 
lowing expression for u^: 

(drai'O + dri^l) F^ 



2Lj^eodzMD + D/l3) 
We stress that, addressing the splitting %p — ipo + ipi, the azimuthal Lorentz force stands as 

F^ ~ [dzl {dr.ijQ + dri>l) - drldz^l] / {^TTrD . (17) 



Expressing ea. ([T5|) for Vr and eg. ([16)) for Vz via the function Q, we get the following integro- 
differential compatibility condition 

dziJi [Jdz drF^] = F^ {dr„^0 + 9,Vi) , (18) 

from which, provided i/'i, i.e., Y{x, u), from eas. dlll) . one can determine the toroidal magnetic 
field J(V',^). 



4. Non-Linear Limit - Let us now analyze the case Y ^ 1, i.e., when the electromagnetic 
reaction within the disk plasma is strong enough to significantly deform the background 
magnetic field. In this context, the vertical and radial equilibria are now described by eq. ljllap 
and 

{D + D//3)Y + \d.,P + [dlY + e,dlY)d^Y = , (19) 

respectively. Clearly, we got a new configuration system for which linear terms in the per- 
turbed magnetic surface Y have been neglected in favor of the quadratic ones. In the same 
limit, the azimuthal equation (1151) rewrites 

2tVrU:K ^F^^ [djdr^i - drld.^i] / {AnrD , (20) 

and eqs. ([T5| . ([T5| rewritten in terms of Y as follows: 

v^^-[vrd^Y]/[^,duY], (21a) 

duY [Jdu d^F^] = F^d^Y , (21b) 

respectively. We now look for a separable solution of the configurational system of the form 

r(^t^ x) = F{u^) sin(ax) , (22a) 

b{u^,x) = d{u^) cos(ax) , (22b) 

F(^t^ x) ^ S,(u^) sin^{ax) + tt{u^) cos{ax) , (22c) 

F^(w2, x) = G{u^) sin'' (ax) . (22d) 

Substituting eqs. (|22[) in the set of equations derived above, we get ordinary differential rela- 
tions in the variable u: from eas. (|llap . p^ . we get the following system 

d^2n + e^d^Q , (23a) 

TT - 2DF/a = , (23b) 

a„2^-f 2(-a2F-|-e2F")9„2i^ = 0, (23c) 

UlFd + aS,+aF{-a'^F + e^F")={) . (23d) 

Using e = e + e, 00. (^0)) and ea. (j21ap give 

Vr = [Gsin^(ax)] / [2ujKto{D + d coii{ax) / P)] , (24a) 

V, = -[aF cotg(ax) w,] / [F'V^] , (24b) 

respectively. Finally, the compatibility equation (|21bp can be recast to give the integro- 
differential relation: 

hF' [jdu G]=FG => [Jdu G] = gF^ , (25) 

where we have used (...)' = d{...)/du and g = const. 

5. Analytical Solution - To derive an exact solution of the system fixed above, we now 
consider the case of large /3- values for the disk plasma, i.e., ez ^ 1. Such analytical solution 
writes 

D = l-ezu'^, i^ = Aexp[-uV2] , A>1, (26a) 

a = V3 , d = Fl{e,Vi) , C = (3 - e.u')F' , (26b) 

TT ^ 2(1 - e^^t2)i^/^/3 , G = -bguF'' . (26c) 



Hence the radial velocity (|24ap takes the explicit form 

1 bgA''ue-''"^/^ sin'' {V3x) 



(27) 



2ujKeo 1 + V3e^Ae-«V2cos(V^x) ' 

where, here and in the following, we neglect the term ezu'^ with respect to terms proportional 
to A. It is remarkable that the total mass density is expressed by the relation e — eo{D+D//3), 
i.e., 

e==eo[l + \/3e^Ae-"'/2cos(V3x)] , (28) 

and e > implies ^/Se^A ^ 1. This feature of the mass density profile (Figure [U Plot (a)) 
ensures that the ring sequence, predicted in [5], is present in this regime too. Since Vr is 




Figure 1: Behavior of e/eo (Plot (a)) and Vz 
A = 100 and e, = 0.005. 



Vz/v* (Plot (b)) {v* = g/iOKeo), setting 



odd in u {i.e., 6 is even in z) the accretion rate of the disk vanishes identically. Despite 
a local non-zero radial matter flux is present in the equilibrium configuration, however the 
compensation between the contributions over and below the equatorial plane, respectively, 
gives a net zero accretion of the disk. 

Finally, the vertical velocity component reads as 



y^SjTz bgA'^e-'"' /^sin''~^(V3a;)cos(^/3a;) 
2ujKeo l + y/3 ezAe-""/^ cos( VSx) 



(29) 



Let now consider the peculiar case 6 = 1 of eq. (|^^ . near the equatorial plane. For a;-\/3 ^ n, 
we deal with a positive or negative vertical velocity (depending on the sign of Ag) having 
a very large amplitude, | Wz |:^ 1 (Figure [H Plot (b)), as far as the condition e^-y/SA ^ 1 
holds {i.e., in correspondence of the mass-density nodes). Such a condition is satisfied only 
in specific layers of the disk configuration since our model is localized at tq but A, g and Cz 
vary with it. 

Finally, the total term describing the magnetic surface, once we account for the first-order 
corrections around tq, reads 



Yt 



Y = X + Acxp[-M^/2] sin(\/3a;) . 



(30) 



The ratios of the magnetic field components to the background magnetic field due to the 
central object, i.e., Bq = Bzo^z, read as 

Bz/BzQ = d^Yt ~ V3Aexp[~u^/2] cos(\/3x) , 

Br/BzQ == -^duYt = y/e^Auexp[~u'^/2]sm{V5x) . 



We stress that the condition to neglect the inertial terms due to the radial and vertical 
velocity (despite of the presence of Vz-peaks) in the system ([TT|) takes the following form: 
Ag/{1 — y/3ezA) ~ C(l). Moreover, in the regime Cz <C 1 (for which our solution is derived), 
the condition to deal with a positive mass density, i.e., \^€zA < 1, ensures the positivity of 
the total pressure too. 

6. Astrophysical Characterization - The configuration scheme we adopted above is appro- 
priate to describe a local radial layer (a small circular crown) of a thin accretion disk since 
all the morphological prescriptions of this astrophysical system are included. The choice of 
the background magnetic surface as a function of tq only, models a dipole magnetic-field of 
the central object, as far as the disk is sufficiently thin and we take uJk oc i/'o- The grav- 
itational field is correctly represented in the thin disk approximation, for the limit of an 
high Keplerian angular velocity to have a gravostatic confinement of the vertical profile of 
the disk, i.e., uj'^z ^ 0{1). The idea that the poloidal (radial and vertical) matter fluxes 
are taken small with respect to the toroidal Keplerian motion is a natural approximation, 
valid in the standard model too. The vertical velocity must be smaller than the azimuthal 
one in a thin disk, because of the background vertical equilibrium, while the radial motion 
is required perturbative since its presence must preserve the steady character of the matter 
distribution. In the standard model the amplitude of Vr is determined by the turbulent vis- 
cosity coefficient, while here it is fixed via the self-consistency of the electron force balance 
equation. Indeed, we deal with a vanishing net accretion, but, from a physical point of view, 
we could speak of a weak constant rate of matter infall, since this fact do not affect our main 
goal about the vertical matter jets. In this ideal plasma scenario, we can reach very high 
self-consistent magnetic fields and we show how they can induce intense axial matter fluxes. 
This picture is commonly believed the basic mechanism of the emission from high energy 
astrophysical sources (such as Gamma Ray Bursts or Active Galactic Nuclei), but its setting 
is still rather questioned. We offer here an alternative explanation, surprisingly relying only 
on the morphology of ideal and highly non-linear MHD structures. 

7. Conclusions - We have shown how the steady state of an accretion disk, described by 
ideal MHD, is compatible with the presence of significant radial and vertical matter fluxes. 
Such matter currents are unable to account for a non-zero accretion rate, but they are relevant 
in view of the instabilities that can be driven by their existence. In particular, we outlined, 
as the main issue of our study, that there are conditions under which strong vertical matter 
fluxes are locally allowed, especially near the equatorial plane. This feature deserves attention 
because dealing with matter jets along the symmetry axis is one of the expected features of 
the disk configurations, but it is hard to be recovered in the standard viscoresistive model. 

** This work was developed within the framework of the CGW Collaboration (www.cgwcollaboration.it).** 
Acknowledgment: NC gratefully acknowledges the CPT - Universite de la Mediterranee Aix- 
Marseille 2 and the financial support from "Sapienza" University of Rome. 



References 

[I] N.I. Shakura, Sov. Astron. 16, 756 (1973). 

[2] G.S. Bisnovatyi-Kogan, R.V.E Lovelace, New. Astron. Rev. 45, 663 (2001). 

[3] E. Velikhov, Sov. Phys. JETP 36, 995 (1959). 

[4] S.A. Balbus, J.F. Hawley, Rev. Mod. Phys. 70, 1 (1998). 

[5] B. Coppi, Plasma Phys. Contrl. Fus. 36, B107 (1994). 

[6] G.I. Ogivile, Mon. Not. RAS 288, 63 (1997). 

[7] B. Coppi, Phys. Plasmas 12, 057302 (2005). 

[8] B. Coppi, F. Rousseau, Astrophys. J. 641, 458 (2006). 

[9] B. Coppi, EuroPhys. L. 82, 19001 (2008). 

[10] B. Coppi, F. Rousseau, in the proceedings of The .34^^ EPS Conference on Plas'ma Physics, Warsaw 
(July 2007), Paper 04.034. 

[II] G. Montani, R. Benini, Mod. Phys. L. A 24, 2667 (2009). 
[12] M. Lattanzi, G. Montani, EuroPhys. L. 89, 39001 (2010). 
[13] D. Lynden-Bell, Mon. Not. RAS 279, 389 (1996). 

[14] H.C. Spruit, Lect. Notes Phys. 794, 233 (2010). 

[15] R. Stehle, H.C. Spruit, Mon. Not. RAS 323, 587 (2001). 

[16] V.C.A. Ferraro, Mon. Not. RAS 97, 458 (1937). 



